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Abstract 



We report on a study of the thermodynamics of the Ginzburg-Landau 

model for two-dimensional type-II superconductors near the transition from 

the normal state to the Abrikosov-lattice state. We couch our analysis in 

terms of the participation-ratio entropy, s(P), which expresses the volume in 

order-parameter-space with a given participation-ratio for the local superfluid 

density. s(P) completely determines the thermodynamics of the system. We 

report on results for s(P) obtained analytically using perturbation expansion 

methods and numerically using Monte Carlo simulation methods and discuss 

the weak first-order phase transition which occurs in this system in terms of 

the properties of s(P). 
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The discovery of high-temperature superconductivity in materials with strong planar 
anisotropy has led to renewed experimental ]1]-|3J] and theoretical |§-|9| interest in the prop- 
erties of two-dimensional and strongly anisotropic three-dimensional type-II superconductors 
at fields near H&iT) where thermal fluctuations are important. Phase transitions of type-II 
superconductors in a magnetic field are unusual because of the Landau level degeneracy of 
Cooper-pair states in a magnetic field flO |. In mean- field-theory continuous phase transi- 



tions occur simultaneously in a number of channels equal to the Landau level degeneracy 
and the low-temperature Abrikosov lattice state has both superconducting and positional 
order. The Cooper-pair Landau level degeneracy also increases the importance of fluctu- 
ations and as a result the phase transition is considerably rounded. It is well established 
experimentally and theoretically that the Abrikosov vortex lattice state in D = 3 dimension 
melts into a vortex liquid state through a weakly first order phase transition, [j3],[4||6|, |iT||L2"|| A 
consensus |p|,p!3|-[T5| has emerged from recent work that a weak first-order phase transition 



with a latent heat which is ~ 2% of the mean-field condensation energy at the depressed 
transition temperature also survives in D = 2, although this view is not universally held 
We have previously |TJj introduced a formulation of the thermodynamics of this system 



in terms of a quantity, s(P), which we refer to here as the participation-ratio entropy and 
which measures the volume in order-parameter space associated with a given participation 
ratio for the local superfluid density. In this Letter we report on evaluations of s(P) based 
on high-temperature and low-temperature expansions and on Monte-Carlo simulations of 
the Ginzburg-Landau (GL) model. 

For fields sufficiently close to H c2 the order parameter is confined to the lowest Cooper- 
pair Landau-level |17]]. The free energy density of the lowest-Landau- level GL (LLL-GL) 



model is 

/[#] = ( a (T) + heH/m*c)\^\ 2 + ^| 4 (1) 

where oc(T) = a (T — T c0 ), T c0 is the zero-field transition temperature, m* parameterizes the 
energy cost of spatial variation of the order parameter and /3 is taken to be independent of 



T. The mean-field theory transition temperature satisfies o.h(T C 2) = ot{T C 2) + heH/m*c = 0. 
The LLL-GL model free energy is F GL = J drf[^f(r)]. For T < T c2 the quadratic term 
in Eq. (|1]) lowers Fql while the quartic term always makes a positive contribution. Our 
approach to the thermodynamics of this system is based on the observation that for a 
given magnitude of the quadratic term in Fgl^] the quartic term is larger when the order 
parameter has a smaller participation ratio ||18|| , 

( fd 2 f\^\ 2 ) 2 

P m = ^ a '\*\ , (2) 

P[^f] is roughly the fraction of the area A of the sample over which the local superfluid density 
(|\1/| 2 ) is spread. (In LLL-GL theory P _1 [\l/] is known as the Abrikosov ratio, For 
a given participation ratio the minimum of Fql/A is —a 2 H P/2(3. Any order parameter in 
the LLL has one zero |]8|JT9| (vortex) for each magnetic flux quantum through the system 
and therefore cannot have the constant value required to obtain P = 1. The mean- field- 
theory order parameter of the LLL-GL model is the LLL order parameter with the maximum 
value of P; P is maximized by placing the vortices on a triangular lattice; Pa = — 
0.862370 • • -. At finite temperatures the properties of the LLL-GL model are determined 
by a competition between the thermal weighting factor (exp(— Fgl\^\I^bT)) which favors 
large values of P[^f] and the distribution function of P[^f] which is peaked at smaller values 
of P as we discuss below. 

We choose to work in the Landau gauge (A = (0, Bx, 0)). The order parameter $?(f) can 



then be expanded in the form, |L4 

*(r0 = (^^) 1/2 E^( 7rl/2 V)' 1/2 exp(^)exp(-(x-H 2 ) 2 /2£ 2 ) (3) 

P k 

where the number of terms in the sum over k is = L x L y /27r£ 2 and £ 2 = h/2eB. We 

define the following intensive thermodynamic variables, which will be used to describe the 
system. A dimensionless average local superfluid density is defined by 

^-k^^l^-w^ (4) 



The participation ratio, defined by Eq. (@) is 

p[C k ] = (^k\Ck\ 2 ) 2 ^ 

J2k!k2k3k 4 Ck 1 Ck 2 Cf C3 Ck 4 5k 1 +k 2 ,k3+k4^ > (ki, k 2 , fc 3 , fc 4 ) 

where 9(A; 1 , k 2 , k 3 , fc 4 ) = (N^LJ Ly) 1 / 2 exp{-§ [E* =1 kf - f (Eii ^) 2 ]}- In terms of A and 
P, the LLL-GL free energy has the following form: 

F GL [C k ] 2/ / \a rr^i i (^o[Cfc]) 2 N 

W = 5 ( s ^)A [C k] + - 1 p^ r )- (6) 



Temperature and field enter only through [fL7] the dimensionless parameter, g 
a H {^ 2 / (3k B Tf' 2 oc (T — T c2 )/{THy/\ 

The partition function for the LLL-GL model is 



J] / dC k dC k exp{-F GL [C k ]/k B T}, (7) 



and can be rewritten as the following form: ||I4] 



Z = J dAodPexp{-iV(A ,P,2 2 )} (8) 

where 

/(A , P, g 2 ) = g 2 (sgn(a H )A + ^) - s(A , P) - Inflat^ 3 //?) (9) 

and 

s(A ,P) = ^lnJI / dC k dC k 5{P-P[C k ])5{A - A [C k \). (10) 
7 V k J 

Because P[Cfc] is invariant under a scale change of |\l/(r)| 2 , it follows that 

S (A ,P) = mAo + s(P). (11) 

The participation-ratio entropy, s(P) = s(l, P), expresses the portion of volume in the phase 
space with a given participation ratio. 

In the thermodynamic limit fluctuations in AofCfc] and P[Cfc] are negligible so that the 
free energy of the system at a given temperature (g) can be obtained by simply minimizing 



f(Ao,P,g 2 ) with respect to Ao and P. It follows that the equilibrium values of Ao and P 
satisfy 

g 2 (sgn(a H ) + ^A ) - ^- = (12) 

and 

*W = -^ A o- (13) 

Eq. ( |T2] ) establishes a functional relationship between Ao and P which originates in the fact 
that the properties of the system do not depend on an and (3 independently but only on 
g oc a H /f3 1/2 . Note that for g » 1 A = while for g « 1 A = 2P. Eq. (0) then 
fixes the equilibrium values of P and Ao and hence the free energy. This equation reflects 
the balance between the rate of increase of volume in order-parameter-space and the rate of 
decrease of condensation energy which fixes the equilibrium participation ratio. Note that 
the right hand side of Eq. flUf ) vanishes for g » 1; the equilibrium value of P in the high 
temperature limit is the value of P where s(P) is maximized, i.e. the most probable value of 
P in order parameter space. We remark a stable equilibrium can occur at any temperature 
at participation-ratio P only if 

s"(P) < 0. (14) 

We now turn to the evaluation of s(P). Since P is positive definite this function is 
defined on the interval (0, Pa) and we can evaluate it analytically near both end points of 
the interval. For P near Pa we can approximate P[C&] by a Taylor series expansion to 
second order around an Abrikosov vortex lattice state. Since the Abrikosov lattice states 
occur at an extremum of P[Cfc] linear terms are absent and by formally diagonalizing the 
resulting quadratic form we find that the volume in order-parameter-space with a given 
value of P is proportional to the surface area of a sphere in ~ 2N^ dimensions with radius 
proportional to (Pa — P) 1 / 2 . It follows that for P near Pa 

s(P) ^ln(P A -P). (15) 



The participation-ratio entropy for small P is most easily evaluated by using the sym- 
metric gauge expansion of the order parameter || in terms of eigenfunctions with definite 
angular momentum, m. In this case order parameters which are confined to an area ~ k/N^A 
centered on the origin, must be expanded in terms of the eigenfunctions with m < k. It 
follows that the order-parameter-space volume with P < k/N^ at A = 1 is given by the 

1 /2 

surface area of a sphere in k dimensions with radius and hence that for P going to zero 

s(P) « —P In P. (16) 
We can also obtain analytic results for the expansion of s(P) about its maximum by 



using the high-temperature expansion of the free-energy of the GL-LLL model pO 



F = N^k B T(\n(-f-) + f 2D (x)) (17) 



where an = 6t x (l — 4x), x = (/3/c^T) /(47r£ 2 o^) and g — (1 — 4x)/v / 4x~- (x is non-negative, 
x = 1/4 at T C 2, x — > for g — > +oo and x — > oo for g — > — oo.) Coefficients in the power 
series expansion of /2£>(x) may be evaluated using a diagrammatic perturbation expansion 
in terms of self-consistent Hartree correlation functions. The six leading coefficients were 



obtained by Ruggeri and Thouless El and the expansion was later extended to eleventh 



order by Brezin, Fujita and Hikami |H] and recently to thirteenth order by us. [pz]] . The 
equilibrium values of A and P may be expressed in terms of /2d(^) by differentiating the 
free energy with respect to an and (3. [^2| ; 



P(x) (l-2s/U*)) s 



(4+(l-4x)/^(x))(l + 4x) 
= - + - x - \x 2 + x 3 - 431.59018759018753 x 4 + 7170.5968205856756 x 5 

-134096.68933891651 x 6 + 2772357.9400791259 x 7 - 62615750.777811569 x 8 

+ 1532019484.1067800 x 9 - 40349691260.478735 x 10 + 1138241888638.9989 x 11 

-34247891779884.099 x 12 + • • • . (18) 

Using Eq. ( |T3| ) and the high temperature expansion of Ao(x) we obtain the following expan- 
sion for the parametric dependence of s'(P) on x: 
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(l-2xf' D (x)) 2 

= -4 x + 8 x 2 - — x 3 + 1020.8000000000000 x 4 - 14012.376334776334 x 5 
3 

+219842.58347280514 x 6 - 3870523.5408992074 x 7 + 75492832.182500102 x 8 
-1615943686.2025482 x 9 + 37680981633.278216 x 10 - 951344015946.40479 x 11 
+25869507417290.368 x 12 + ■ • • . (19) 

Eq. (|18D can be inverted to expand x in terms of P = P — 1/2. Inserting this series in 
Eq. ( |T9"|) gives s'(P) as a series in P. Integrating this series and determining s(P = 1/2) by 
direct evaluation we find 

160 ~ 1096 ~ 

S(P) = l + lnvr + 0P-4P 2 P 3 + P 4 - 802.2613949013954 P 5 

K ' 9 45 

+ 16926.259949713371 P 6 - 499062.26189198034 P 7 + 17336002.792339820 P 8 
-688621672.74816198 P 9 + 30506759410.693500 P 10 - 1483164272081.5901 P 11 
+78249334359919.989 P 12 - 4439326835886646.1 P 13 + ■ • • . (20) 

Directly evaluating the low-temperature expansion of the LLL-GL model has proven to 
be an arduous task fl23|-f2"5H and has led to discordant results. In contrast, using Eq. ([15]) 



we can obtain a result for the the leading low-temperature correction to the mean-field free 
energy of the LLL-GL model in a very simple way. We find immediately, in agreement with 
Ref. jpS], that for large x (g « 1 ) 



, / x (l-4x) 2 P A , 

f 2D (x) = -± ^ + lnx. 21 

4x 

The first term here is the mean-field-theory energy and the correction comes from the low- 
temperature equilibrium participation-ratio entropy. 

It is interesting to observe that equilibrium values of P vary through a relatively narrow 
range between the largest possible value for P[^] (Pa) and the most probable value for P[^] 
(1/2) from low temperature to high temperature limits. Useful approximate expressions for 
the free energy, the magnetization, and the specific heat of the GL-LLL system have been 
proposed by Tesanovic and collaborators [p6jp7| motivated by this property. 



The high-temperature expansion of s(P) around P = 1/2 can be extrapolated with 
the use of Pade approximants for s'(P). The Pade approximants were chosen to satisfy 
s'(P) = 1/(P — Pa) for P — > Pa- and s(P) was obtained by integration. Results are shown 
in Fig. [I]. Poles appear in the approximants for P ~ —0.2 and the extrapolation to negative 
values of P is not very successful. Nevertheless, we believe that s(P) is a smooth function 
over the entire interval (0, Pa)- This expectation is consistent with numerical results for s(P) 
obtained by Monte Carlo methods which are also shown in Fig. |l[ The Monte Carlo results 
were extracted from distribution functions for participation-ratio values, A\(P), calculated 
using exp(— A^(AP + A )) as the sampling function. Since A\(P) oc exp(N^(s(P) — AP)), 
extrema of the distribution occur where s'(P) = A. By performing calculations at a series of A 
values we were able to map out the function s'(P); the results shown in Fig. p] were obtained 
by numerical quadrature from the Monte-Carlo results for s'(P). The overall agreement 
between the analytic and numerical results for the participation-ratio entropy is excellent. 
The inset in Fig. [I] shows Monte-Carlo results for s'(P) in the narrow range of equilibrium 
participation-ratio values which occur near the first-order melting transition. We see that in 
the Monte-Carlo simulations (but not in the analytic results) s"(P) > for P in the interval 
(0.832, 0.837). Equilibrium P values cannot occur in this interval and must therefore have a 
discontinuity in their temperature dependence. It is this property of the participation-ratio 
entropy which leads to the weakly first order phase transition in the LLL-GL model. 

Our description of the thermodynamics of the LLL-GL model is summarized in Fig. |2| 
in terms of three contour plots for f(A ,P,g 2 ) + ln(|af/|7r£ 2 / ' 0) . At each g the equilibrium 
(A ,P) minimizes f(A ,P,g 2 ). The top panel is for a temperature at which the system is 
in the vortex liquid state, the middle panel is for a temperature close to the phase transition 
and the bottom panel is for a temperature at which the system is in the vortex lattice state. 
The trend to decreasing participation ratios at higher temperatures is driven by the increase 
in the relative importance of the entropy. The first order phase transition occurs because 
of the occurrence of an interval over which s"(P) is positive. When the optimal P values 
are close to this interval two local minima appear in the free energy contours and the global 
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minimum switches from the local minimum at larger P to the local minimum at smaller P as 
the temperature increases. In the high temperature limit of the model, thermal fluctuations 
are Gaussian and P approaches 1/2. 

This work was supported by the Midwest Superconductivity Consortium through D.O.E. 
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FIGURES 



FIG. 1. Analytic and numerical data for s(P). The inset shows Monte Carlo results for s'(P). 

FIG. 2. Contour plots for f(Ao,P,g 2 ) + ln(|a/f|7r£ 2 //3) at g = — 5.5(top panel, vortex liquid 
state), g = —6.6 (middle panel, at phase transition ) and g = —7.1 (bottom panel, vortex solid 
state) . 
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